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ABSTRACT 

The young, low-mass, triple system NTTS 155808-2219 (ScoPMS 20) was 
previously identified as a ~17-day period single-lined spectroscopic binary with 
a tertiary component at 0.21 arcseconds. Using high-resolution infrared spectra, 
acquired with NIRSPEC on Keck II, both with and without adaptive optics, we 
measured radial velocities of all three components. Reanalysis of the single-lined 
visible light observations, made from 1987 to 1993, also yielded radial veloc- 
ity detections of the three stars. Combining visible light and infrared data to 
compute the orbital solution produces orbital parameters consistent with the 
single-lined solution and a mass ratio of q = 0.78 ± 0.01 for the SB. We dis- 
cuss the consistency between our results and previously published data on this 
system, our radial- velocity analysis with both observed and synthetic templates, 
and the possibility that this system is eclipsing, providing a potential method for 
the determination of the stars' absolute masses. Over the ~20 year baseline of 
our observations, we have measured the acceleration of the SB's center-of-mass 
in its orbit with the tertiary. Long-term, adaptive optics imaging of the tertiary 
will eventually yield dynamical data useful for component mass estimates. 
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Introduction 



To understand the star formation process, reliable observations of fundamental stellar 
properties are needed so theoretical models can be tested. From on e star forming region 



(SFR) to the next the bin ary fraction can vary ( Mathie u et al.ll2000l ) and may depend on 



the density of the region ( jSimonl 119971 ; iBeck et al. 



20031 ) and spectral type of the primary 



( jLadal 120061 ). underscoring the importance of a complete binary census. Characterization of 
multiples' frequency, component separation distribution, and mass ratio distribution for a 
given SFR provides clues to the broad star forming properties (angular momentum, density, 
turbulence, etc.) of the parent molecular cloud. Small separation spectroscopic binaries 
with periods sufficiently short to enable the measurement of the individual stellar velocities, 
and thus of the system's mass ratio, are p otential targets for the dynamical determination of 



indiv idual component stellar masses (e.g. JSteffen et al.ll200ll ; iPrato et al.ll2002al ; iBoden et al. 



20051 ). Knowledge of absolute masses and observable properties, such as effective temperature 
and l uminosity, plays a k e y role in the improvem ent of pre-main sequence evolutionary models 



Baraffe et al 



1998 



Palla fc Stahlerll200lf ). This work is part of our effort to measure a 



(.e.g., 

substantive sample of young star mass ratios and masses with the ultimate goals of improving 
our understanding of tight binary formation and of tying models of young star evolution to 
concrete dynamical data. NTTS 155808—2219 is a particularly important target given its 
high-order multiplicity and its low-mass components, a regime plagued by discrepancies in 
theoretical calculations. 

Spectroscopic binaries (SBs) with low mass-ratios are usually identified as single-lined 
spectroscopic systems when observed in visible light because the large difference in flux 
between the primary and secondary at short wavelengths prevents detection of the secondary 
component (Mazeh et al. 2002). In the Raleigh- Jeans regime of the stars' spectral energy 
distribution, however, flux scales much less steeply as a function of mass. Thus, by observing 
single-lined spectroscopic binaries with infrared (IR) spectroscopy we are able to improve our 
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Observations of Einstein x-ray sources in Upper Scorpius by IWalter et al.l (I1994J . here- 
after W94) identified the spectral type M3 NTTS 155808-2219 system as an SB with a 
tertiary component (SB+T). Figure 11 of W94 presents a visible light spectrum with three 
distinct Li A6707 A lines. Based o n this evi d ence for multiple components, they estimated 
an SB mass ratio of 0.71 ± 0.05. iMathieul (119941 ) provided initial single-lined orbital pa- 



rameters for the NTTS 155808— 2219SB, revealing a 16.925 day period, -5.0 kms -1 center- 
of-mass veloc i ty (7 ), consistent with membership in Upper Sco, and an eccentricity of 0.10. 
Kohler et al.l (120001 ) later determined an SB+T separation of 0.193", which is consistent with 
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the value of 0.21" determined by I Woitas et al.l ( 2001 ). Equivalent widths of the 6707A Li 



line (~0.6A) and Ha (~ -3A) suggest a young age and a low accretion rate for the system. 
The lack of evidence for circumstellar dust, based on 3.4-22 fim photometry from the Wide- 
field Infrared Survey Explorer (WISE), is consistent with a weak-lined T Tauri star (WTTS) 
classification. Observed properties of NTTS 155808—2219 are presented in Table [TJ 

This paper describes the results of combining high-resolution IR spectroscopy with re- 
processed visible light data to determine the component radial velocities of NTTS 155808—2219 
over a nearly 20 year timespan. In §2 we briefly describe our observations and data reduction. 
Our cross-correlation analysis appears in §3, and the NTTS 155808—2219 orbital solution is 
presented in §4. Section 5 contains our discussion of the system age and component masses, 
the possibility of an eclipse, given the large SB orbital inclination, and initial estimates for 
the tertiary orbital characteristics. A summary of our findings appears in §6. 



2. Observations and Data Reduction 

2.1. Visible Light Spectroscopy 

Spectroscopic observations of NTTS 155808—2219 were collected from 1987 June to 
1993 April on the 37 epochs listed in Table |2] (W94). Two nearly identical echelle spec- 
trographs, with photon-counting intensified Reticon detectors, were used on the Multiple 
Mirror Telescope on Mount Hopkins, Arizona (prior to its conversion to a monolithic mir- 
ror), and on the 1.5-m Tillinghast reflector at the F. L. Whipple Observatory, also atop 
Mount Hopkins. A single echelle order, 45 A wide, was recorded at a central wavelength of 
5188.5 A, containing the Mg I b triplet. The resolving power provided by this setup was 
R= A/AA ~ 35,000. The nominal signal-to-noise ratios of the 37 spectra obtained range 
from ~6 to 19 per resolution element of 8.5 kms -1 . Standard flatfield frames were obtained 
each night, and the pixel-to- wavelength mapping was b ased on exposures of a thorium- argon 



lamp taken before and after each science exposure (see lLatham et al.ll2002l ). 



2.2. Near-Infrared Spectroscopy 

Observations were made in the IR between 2000 June and 2007 April with the Keck II 
10-m telescope on Mauna Kea. The UT dates of observation are listed in Table [3J Nine H- 
band orders, with a central wavelength of ~1.555 fim in the middle order, were obtained using 



the f a cility near-infrared, cross-dispersed, cryogenic spectrograph NIRSPEC (IMcLean et al. 



19981 . 120001 ) . NIRSPEC employs a 1024 x 1024 ALADDIN InSb array detector. Source 
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acquisition was accomplished with the slit viewing camera, SCAM, which utilizes a 256 x 
256 HgCdTe detector. NTTS 155808-2219 has a 2MASS if-band magnitude of 9.01 (Table 
1); integration times for individual frames were 300 seconds with the NIRSPEC-5 filter. 
We nodded the telescope to move the star between two positions on the slit, separated by 
~60 pixels, to allow for background subtraction by differencing sequential spectra. 

Four of the seven observations used the 24" x 0.288" slit and contain light from the 
three SB+T components. The Keck Adaptive Optics (AO) system effectively reimages the 
optics, reducing the plate scale by a factor of ~11, producing slit dimensions of 2.3" x 
0.027". Three IR observations utilized AO in front of NIRSPEC. On all three AO nights we 
obtained isolated spectra of the SB. On two of these occasions we utilized SCAM imaging 
to align the slit to also include the tertiary. All spectra have resolution R w 30,000. The SB 
and SB+T spectra from the seven epochs of observation are shown in Figure [T] for a single 
H-band order. 

All spectroscopic reductions were made by using the REDSPEC package, software writ- 
ten at UCLA by S. Kim, L. Prato, and I. McLean specifically f or the analysis of NIRSPEC 
dat£0, following the procedures outlined by lPrato et al.l (J2002bJ). We used the central order, 



49, found at our setting. This order has three advantages: (1) it is rich in both atomic 
and molecular lines and is therefore suitable for identifying spectra of both warm and cool 
stars, (2) the OH night sky emission lines across order 49 are numerous and well-distributed, 
yielding accurate dispersion solutions, (3) this order has the advantage of lacking prominent 
telluric absorption lines. Consequently, we did not have to divide our target spectra by 
telluric star spectra. Because the AO reimaged slit width is so narrow, it was not possible to 
use OH night sky emission for wavelength calibration on the AO nights. On these occasions 
we took comparison lamp exposures at the beginning and/or end of the night to determine 
the dispersion solution and zero point. 



3. Analysis 

3.1. Visible Light Radial Velocities 

Based on the expectation that the secondary would be too weak to detect in our spectra, 
and that only the primary and tertiary components would be visible, we initially derived 
radial velocities for only these two components of the NTTS 155808—2219 sy stem using 
TODCOR, a two-dimensional cross-correlation technique (IZucker fc Mazehlll994l ). The same 



http://www2.keck.hawaii.edu/inst/nirspec/redspec/index.html 
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synthetic te mplate was used for the primary and tertiary stars, based on Kurucz model 
atmospheres (ILatham et al.ll2002l ). with a temperature of 3750 K selected in accordance with 
the late spectral type of the unresolved system, M3, and zero rotational broadening, vsim = 
kms -1 . The large scatter in the tertiary velocities suggested possible contamination either 
from moonlight in a few exposures, or more likely by th e secondary compon ent. A reanalysis 
using an extension of TODCOR to three dimensions ( jZucker et al.lll995l ) showed that the 
lines of the secondary are indeed visible in most of our visible light spectra, and radial 
velocities were derived using a template for the secondary identical to that of the primary 
and tertiary The moonlight contaminatio n was removed from seven of the exposures using 
a four- dimensional version of TODCOR ( {Torres et al.l 120071 ) and the addition of a fourth 
template appropriate for the Sun. The radial velocities are listed in Table |2l and have 
typical uncertainties of 1.7 kms -1 for the primary, 4.4 kins" 1 for the secondary, and about 
3.5 kms -1 for the tertiary. The stability of the zero-point of the CfA velocity system was 
monitored by means of exposures of the dusk and dawn sky, and small systematic run-to- 
run corrections were applied in the manner described by iLathaml ( 119921 ). already included 
in the velocities described above. In the case that either the secondary or the tertiary radial 
velocity was impossible to determine, only the primary star velocity was used. This resulted 
in 37 primary, 33 secondary, and 33 tertiary velocities. The light ratios among the three 
stars determined from our spectra are a = = 0.30 ± 0.03 and (3 = izfl\ = 0.53 ± 0.04 
at the mean wavelength of our observations. The tertiary companion is thus brighter than 
the SB secondary. Comparis on of the visible light spectra to GL 752A, a main sequence 
star with vsim < 2.5 kms -1 ( [Browning et al.ll2010l ). implies a low, < 10 kms -1 , rotational 
velocity for all components of NTTS 155808-2219. 



3.2. Infrared Radial Velocities 

Stellar radial velocities were also measured for the IR spect ra using a simila r multi- 



dimensional cross-correlation analysis following the approach of IZucker fc Mazehl (11994 ). 
To optimize the correlation, and thus obtain the most precise radial velocities, the best 
matching synthetic or observed templates must be used for cross-correlation with the target 
observation. In our analysis of the IR data we tested three sets of templates: (1) the 
NTTS 155808—2219 tertiary spectrum, (2) synthetic spectra, and (3) observed spectra of 
radial velocity standards. 

The tertiary spectra isolated in the AO observations on UT 2000 Jun 11 and 2003 Apr 
20 provided an excellent template because this bound companion likely shares the same 
age, metallicity, and surface gravity with the close SB components. Assuming a distance 
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of 145 pc to the region, the ~0'/21 angular separation of the tertiary implies a projected 
physical separation of ~30 AU. This ability to angularly resolve the subarcsecond, visual 
companion illustrates a unique advantage of IR observations. Cross-correlation of the AO 
tertiary spectra against on e another revealed no significant chan ge in radial velocity, within 
the ~lkms _1 uncertainty (IPrato et al.ll2002bl ; ISteffen et al.ll200ll ). Thus the tertiary spectra 
were combined to use as a cross-correlation template with the highest possible SNR. 

The observed template that best matches the NTTS 155808—2219 tertiary IR spectrum 
is the m ain-sequence star GL 752 A. In a comparison of a sequence of high- resolution, near-IR 
spectra, iPratol (2007 ) found a spectral type of M2 for GL 752A, similar to the M2.5 type of 
Browning et al.l ( 120101 ) based on visible light observations; SIMBAD lists a type of M3. We 
adopt M2 throughout our analysis. Figure [2] shows the GL 752A spectrum together with the 
best matching synthetic template and the average of the two NTTS 155808—2219 tertiary 
spectra taken on AO nights. 

Our synthetic spectral library was gen erated from atmospher ic structures calculated 
with the PHOENIX model atmosphere code ( lHauschildt et al.lll999l ) and kindly provided by 
T. Ba r man. Synthetic templates were compared to observed templates by cross-correlation. 
Pratol (120071 ) determined the spectral types of a subset of ou r observed temp l ates; the corre- 
sponding effective temperatures (T e ff) were estimated from iLuhman et al.l ( 20031 ) . For the 
observed templates, the highest correlations were obtained with syntheti c templates of typ i- 
cally ~500 K hotter T e ff, a known characteristic of model atmospheres (I Jones et al.l 120051 ). 
Because we derive the radial velocity from the cross-correlation, this T e ff discrepancy doesn't 
impact the radial velocities measured with the synthetic templates. Four synthetic templates 
were ultimately used, T e ff = 3645, 3883, 4209, and 48 31 K, corresponding to 1 Gyr old mod- 
els of 0.5, 0.6, 0.7, and 0.8 M stars with logg of ~4.7 faaraffe et al.lll998h . The T eff = 4209 
K synthetic spectrum yielded the highest correlation with the NTTS 155808—2219 tertiary. 

The GL 752A spectrum was c onvolv e d wit h a series of kernels to mimic a range of 
rotational velocities, as described in IPratol ( 120071 . section 4.1). The highest correlation with 
the NTTS 155808—2219 tertiary was obtained for usini = 22 ± 1 km/s. However, radial 
velocities derived with the GL 752A template rotated to other values of vsim vary by a few 
km/s. Cross- correlation of the T e ff — 4209 K synthetic template against the tertiary resulted 
in correlation coefficients that are ~5% lower than for the observed GL 752A template, but 
yielded radial velocities consistent over a wide range of rotational velocities. The highest 
correlation was obtained with the T e ff — 4209 K spectrum rotated to 24 km/s. 

All three types of templates were used to analyze the 2003 April 20 and 2004 May 
24 SB spectra, for which the component radial velocity separation was maximum and the 
SB isolated. Using the averaged tertiary spectrum as a template for both the primary and 
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secondary resulted in the highest correlation. GL 752A as a template for the primary and 
secondary produced nearly the same correlation as the tertiary as the template. However, 
for a range in vsim's, the most consistent velocities were again found using the synthetic 
templates. When the various templates were applied to the SB+T spectra, only the synthetic 
templates produce reliable radial velocities, similar flux ratios to the SB-only analysis, and 
high correlations. 

Testing different synthetic templates to find the maximum correlation confirmed that 
the T e ff = 4209 K spectrum produced the best match to the primary and tertiary, while 
the secondary was best fit by the T e ff = 3883 K synthetic spectrum. Vsim values of 20 
± 1 kms -1 , 24 ± 4 kms -1 , and 24 ± 3 kms -1 for the primary, secondary, and tertiary, 
respectively, yielded the highest correlations. Vsim uncertainties are the standard deviation 
of the highest correlated templates for all IR epochs. Our visible light cross- correlation 
analysis yielded vsim values of < 10 kms -1 . Because the visible light synthetic templates 
used in this analysis were created with inherent rotation, rather than through convolution 
with a kernel function after the spectra were synthesized, we suspect that the lower vsim 
values are likely more accurate. We therefore adopt vsim=10 km/s for all three components 
of NTTS 155808-2219. 

The radial velocities determined with the synthetic template spectra are listed in Tabled 
and have typical uncertainties of 1.9 kms -1 for the primary, 2.6 kms -1 for the secondary, 
and 1.8 kms -1 for the tertiary. The uncertainty in each radial velocity was determined as a 
function of the FWHM of the correlation peak, and the ratio of the cor relation peak heigh t 



to the amplitude of the antisymmetric noise, following the procedure of iKurtz et al.l (119921 ). 
Ideally we would use a three dimensional cross-correlation for all SB+T spectra but the large 
flux ratios and similar spectral types prevented convergence on a solution. In these cases we 
used a two dimensional cross-correlation to fit the primary and tertiary components, then 
identified the secondary radial velocity as a residual peak in the cross-correlation function. 

The three pure SB spectra observed with NIRSPEC+AO yielded an H-band light ratio 
of the secondary /primary, a = 0.55 ± 0.10. For the four SB+T spectra, a was held at 0.55 
during cross-correlation. The 2002 July 18 and 2004 January 27 SB+T spectra were used 
to determine the tertiary/primary light ratio, (3 = 0.81 ± 0.09, and fo r cross-correla t ion o f 



the four SB+T spectra (3 was held at 0.81. H-band magnitudes from IWoitas et al.l (120011 ) 
produce the photometric H-band brightness ratio T/SB = 0.58 ± 0.03, consistent with our 
T/SB H-band value of 0.52 ± 0.13. 
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SB Orbital Parameters 



Prelimin ary estimates fo r the center-of-mass velocity and mass ratio were derived using 
the method of I Wilson! (Il94lf ). Figure [3] shows the primary velocity vs. the secondary velocity, 
for the epochs from which both were determined for the double-lined spectra, along with the 
weighted linear fits for the visible light and IR velocities, separately. From the linear fit, the 
negative of the slope gives the mass ratio, q. A very similar value, to within la, was obtained 
for the visible light and IR data sets. The center-of-mass velocity 7 is given by the y-intercept 
divided by (1+q). Although the full orbital solution provides more precise measurements of 
these quantities by adding phase constraints, the vertical offset between the two data sets 
produces a potentially significant measurement of the center-of-mass acceleration. For the 
13.18 years between the mean epochs of observation d7/dt = -0.30 ± 0.07 (kms _1 yr _1 ). 



Table H] lists the orbital parameters derived by iMathieul (ll994f ) and W94, along with 
the solutions determined using our visible light and IR radial velocities. The orbital fit to 
the radial velocities of both SB components is ac complished with a standard least-squares 
analysis using the Levenberg-Marquardt method ( jPress et al.lll992|). Initial guesse s for the 
solution were found by using an amoeba search routine also from lPress et al.l (119921 ) . Refine- 
ments to the initial guess orbital elements were made to ensure that the best-fit, in the x 2 
sense, was consistent with the single lined solution of IMathieul ( 1l994l ). We determined the 
period, P, the projected semi-major axes of the components, c^sim and a 2 sim, eccentricity, 
e, periastron angle, u, the time of periastron passage, T, the center-of-mass velocity, 7, and 
the center-of-mass acceleration, d7/dt. From the fit to the phased radial velocity curves 
(Figure H]) we can measure the semi-amplitude of both components, K\ and K 2 , and derive 
the mass ratio, q = Ki/K 2 = c^sim / c^sim, and the minimum masses, Misin 3 i and M 2 sin 3 i 
Uncertainties calculated by the Levenberg-Marquardt method are the formal errors from the 
nonlinear least-squares fitting. The orbital elements derived using all the data are dominated 
by the more numerous visible light radial velocities and are consistent with the individual 
IR and visible light orbital solutions. 

On the basis of the combined visible and IR data we find that the SB in NTTS 155808—2219 
has an orbital period P = 16.9243 ± 0.0002 days, center-of-mass acceleration d7/dt = -0.25 
± 0.04 (kms _1 yr~ 1 ), and eccentricity e = 0.113 ± 0.008. Our center-of-mass velocity, 7 = 
-8.06 ± 0.29 kms -1 is cons istent with the valu e of 7 = -6.28 ± 3.04 kms" 1 for Upper Scor- 
pius M dwarfs surveyed by iDahm et al.l (120121 ) .The combined orbital solution, as a function 
of phase, is plotted with the measured radial velocities in Figure HJ In order to depict the 
combined orbital solution as a function of phase, 7 has been subtracted from each of the 
measured radial velocities. The mass ratio, derived from the ratio of projected semi-major 
axes, is q = 0.78 ± 0.01, consistent with the preliminary analysis based on Figure [31 The 
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drift in the center-of-mass velocity of the binary is illustrated in Figure |5l where we have 
subtracted the orbital motion from the individual velocities, and represented the weighted 
average of the residuals as a function of time. The dashed line corresponds to the change in 
7 from the joint orbital fit. 



Discussion 



5.1. Component Masses and System Age 



The template analysis in §3.2 revealed that the observed GL 752A template generally 
produced the highest correlation for all three components and is best matched to the primary, 
which we classify as an M2±l. Cross-correlation using GL 436 (M3) as the secondary 
template resulted in a higher correlation than GL 752A for the 2003 April 20 SB spectrum. 
The secondary was also fit better by a lower temperature synthetic template, resulting in 
our classification of the secondary as an M3±l. Although the GL 752A and T e ff = 4209 K 
templates provided the highest correlation coefficients for the tertiary, the T e ff — 3883 K 
and GL 436 templates were similar to within a few percent. Thus, we classify the tertiary 
at an intermediate spectral type of M2.5±l. Although the synthetic templates yielded more 
consistent RVs, because of the known discrepancy in T e ff, i-e. the T e ff values of synthetic 
spectra are anomalously hot, the observed templates provided a more accurate method for 
estimating the NTTS 155808—2219 spectral types and hence T e jjs. 

To place NTTS 155808—2219 on an H-R diagram, we use t he spectral typ e s M2, M3, 
and M2.5 (§3.2) and the temperature scale for young stars from iLuhman et al.l ( 120031 ): we 
found T e ff = 3560 ± 145, 3415 ± 145, and 3488 ± 145 K for the p rimary, secondary an d 



dc Zccuw et al. 1999). 



tertiary, respectively. We assume a distance to Sco-Cen of 145 ± 15 pc 

Estim ates of A u in the literature are 0.3 (W94) and 0.5 magnitudes (IPreibisch fc Zinnecker 



Prato et al. 



1999h . Using 2MASS magnitudes and the color-color method described in §3.1 of 
fl2003h yields A. u = 0.47. We adopt K u = 0.5±0.2 in our analysis. Finally, using data from 
2MASS and the H-band flux ratios determined by cross-correlation (§3.2), we calculated the 
luminosity of the primary, log(L!/L ) = -0.56±0.04, secondary, logfL^/ L^) = -0.84 ± 0.04 , 
and tertiary, log(L 3 /L Q ) = -0.66 ± 0.04 following the approach in §3.3 of iPrato et al.l (120031 ) 



and using d=145 pc. The NTTS 155808—2219 system lacks any detectable near-IR excess 
(i.e. Tk=0; Prato et al. 2003). Table [5] lists these derived properties. 

The positio ns of the NTT S 1558 08—2219 components on the pre-main sequence evolu- 
tionary tracks of iBaraffe et al.l (119981 ) are shown in Figure EJ From this plot we are able to 
estimate the masses of the primary, Mi = 0.50 ± 0.10 M , secondary, M2= 0.38 ± 0.10 M , 
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and tertiary, M3 = 0.44 ± 0.10 M . These masses give a secondary-to-primary mass ratio 
M 2 /M! = 0.76 ± 0.14, consistent with our orbital solution, a tertiary-to-primary mass ratio 
Ma/Mi = 0.88 ± 0.14, and a tertiar v -to-SB mass ratio M T /M SB = 0.50 ± 0.1 7. Using the 



K-ban d flux ratio from lKohler et al.l ( 120001 ) as a proxy for the mass ratio, as in lKraus et al- 



to 0081 ). results in a T/SB mass ratio of 0.57 ± 0.10, which is in agreement with our estimate. 
We comb ined these mass es timates and our orbital solution (Table @J and followed the ap- 
proach of iPrato et al.l (j200lh to estimate the orbital eccentricity of the system, > 61 degrees 
(Figure E]). 

From the H-R diagram we estimate an age of 3.3 If 3 Myr. All three components of 
NTTS 155808—2219 fall very nearly on the same i sochrone, well within la of e ach other. 
Without adjusting for the multiplicity of the system, iPreibisch &: Zinneckerl ( 119991 ) estimated 
a system luminosity of log(L/L ) = -0.398. Using the initial visible light spectral type of 
M3 fr om W94, and the corresponding temperature T e ff = 3451 K, IPreibisch fc Zinnecker 
( 119991 ) derived an age estimate of 0.5 Myr. Taking the effects of multiplicity into account, 
they derived a luminosity error that increased their age estimate to 1-3 Myr, consistent with 
the age range that we have determined. The upper limit of our age estimate is also consistent 
with the 5-6 Myr ag e determined from the B star main sequence turnoff in Upper Scorpius 
rtde Geus et aDll989l ). 



We use the primary star luminosity and T e ff to derive the radius, Ri 



1 OO+0.19 -r 

-i.oo_ 017 rt , 



which is shown in Figure [8J Curves of Rsini are overplotted for values of v sim of 10 and 
20 km/s. The resulting Rsim curve for fsini = 20 km/s does not intersect with the derived 
radius of 1.38 R , suggesting that the lower value for vsim is more accurate. 



5.2. Possibility of SB Eclipse 

To estimate the probability of at least a partial eclipse, we assume that the primary and 
secondary stars have equal radii, 1.38±q;];7 R , and use the value determined in our orbital 
fit for the semi-major axis, asini = 25.12 ± 0.25 R (Table HJ). For the plausible range of 
stellar radii we expect a partial eclipse for inclinations > 83°, thus for an orbital inclination 
of >61 degrees (Figure 7) we derive an eclipse probability of ~24%. 

Dedicated photometric monitoring of NTTS 155808—2219 could reveal eclipses, and 
thus a value for the orbital inclination, providing absolute masses. The absence of eclipses 
would also be useful in limiting the orbital inclination of the system and thus our mass 
estimates. Bright ness variations aris ing from the spotted primary star's rotation, AV = 



0.164 magnitudes (lAdams et al.lll998l ). might complicate photometric observations of partial 



eclipses, although spot- and eclipse-induced variability have distinct signatures. Another 
difficulty in measuring an eclipse is the nearly integer period of the SB orbit, ~17 days, 
resulting in limited phase coverage for a single observing season. 



5.3. Tertiary Orbital Information 



For the tertiary we h ave two astrometric measurements from the literature, ( IKohler et al. 



20001 ; IWoitas et al.ll200ll ). and another pair determined from our NIRSPEC+AO observa- 
tions. Although the saved SCAM images taken on our AO nights were overexposed and 
saturated, the finite width of the slit and the separated SB and tertiary continua allow us to 
estimate the position angle§ and separation for two epochs. Table [6]lists the astrometric data 
for the tertiary, which show no significant change in the separation and only a small change 
in the position angle over the nearly 9 years of observation. Assuming a distance to Sco-Cen 
of 145 pc, the average on sky separation of 0'.'21 results in a projected separation of ~30 AU. 
For this separation, the masses derived in §5.1, and assuming an edge-on, circular orbit, the 
space velocity of the tertiary relative to the SB center-of-mass would be ~5.5 kms -1 and 
the minimum orbital period would be ~143 years. Averaging our IR tertiary velocities and 
computing the velocity separation between the tertiary and the SB center-of-mass produces 
ARV = -5.73 ± 0.44 kms -1 at the median epoch of our IR data. This velocity separation 
suggests that the tertiary orbit is either circular with a high in clination (implying near copla- 



narity with the SB) or eccentric with any possible inclination; iMonin et al.l (120071 ) show that 
systems with higher multiplicity are more likely to be nonaligned. Although we determined 
d7/dt = -0.25 ± 0.04 (kms -1 yr -1 ) from the SB orbital solution, no significant trend is found 
in the tertiary radial velocities (df^/dt = 0.07 ± 0.09 kms -1 yr -1 ). One explanation for a 
missing trend might be that the tertiary is itself a binary, and thus more massive than we 
expect from a M2.5 alone. Yet, we see no evidence in the tertiary spectrum for another 
star. An alternative explanation for the d7/dt that we have measured is a mismatch of the 
instrumental zero points. However, this is unlikely since both the visible and IR spectra were 
reduced using synthetic spectral templates and the Earth's atmosphere as the radial velocity 
standard frame of reference. Besides monitoring the SB for eclipses, other future work will 
include occasional astrometric imaging and NIRSPEC+AO spectroscopy to separate the SB 
and tertiary and monitor their radial velocities. 



See: http://www2.keck.hawaii.edu/inst/tools/skypa/ 
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6. Summary 



NTTS 155808-2219 (ScoPMS 20) is a WTTS spectroscopic binary (SB) with a tertiary 
component at 0'/21. By combining radial velocities derived from visible light and infrared 
spectra, acquired between 1987 June and 2007 April, we have determined the orbital param- 
eters of the SB. Our orbital solution results in a 16.9243 ± 0.0002 day period, center-of-mass 
acceleration d7/dt = -0.25 ± 0.04 (kms _1 yr _1 ), eccentricity e = 0.113 ± 0.008, and mass 
ratio q = 0.78 ± 0.01. 

Through the use of adaptive optics in the near-IR we were able to disentangle the SB 
and tertiary spectra. Synthetic templates produced consistent and precise radial velocities, 



both in visible light and in the IR. Observed IR templates (e.g., iPrato et al.ll2002bf ) provided 
estimates for the component spectral types of M2, M3, and M2.5 for the primary, secondary, 
and tertiary, respectively. Based on t hese spectral types, effective temperatures of 3560, 
3415, and 3488 K were estimated from lLuhman et al.l ( 120031 ). 



We found H-band flux ratios of a = 0.55 ± 0.10 and (5 = 0.81 ± 0.09, which we used to 
convert near-IR magnitudes t o component luminosities. Plotting the stars on the theoretical 
tracks of iBaraffe et al.l ( 19981 ) gave an age for the coeval system of 3.3^i3 Myr. Comparison 
of the minimum mass from the orbital solution to the H-R diagram mass estimates limits the 
SB orbital inclination to > 61°. Combining T e ff and L, we estimated the primary radius to 
be 1.38lo'}7 R©- Assuming a similar secondary radius, and taking into account the derived 
semi-major axis and limits on the SB orbital inclination, we predict a ~24% chance of a 
partial eclipse of the SB. High cadence imaging of the system would reveal an SB eclipse 
which could help determine the inclination and thus yield component mass estimates. 
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Table 1. Observed Properties of NTTS 155808-2219 



Property 


Value 


R.A. a (J2000) 


16:01:05.2 


Dec. a (J2000) 


-22:27:31.2 


V b (mag) 


13.74 ± 0.03 


J a (mag) 


9.74 ± 0.03 


H a (mag) 


9.01 ± 0.02 


K a (mag) 


8.75 ± 0.02 


L /c (mag) 


8.46 ± 0.09 


M c (mag) 


7.85 ± 0.35 


Wl d (mag) 


8.64 ± 0.02 


W2 d (mag) 


8.53 ± 0.02 


W3 d (mag) 


8.43 ± 0.03 


W4 d (mag) 


8.21 ± 0.28 


Ha EW b ' e (A) 


-3.32 ± 0.72 


Li EW b (A) 


0.55 ± 10% 



a 2MASS All-Sky Point 
Source Catalog 



Walter et all f ll994h 



Jensen fc Mathieul (119971 ) 



d WISE All-Sky Data Release 



Mathieul (11994T ) 



Table 2. Visible Light Spectroscopic Observations and Radial Velocities 



UT Date 


HJD 


7 U 




7?o 


rj 


u tert 


SB 




(2 400 000+1 


( lr m q~ 


( lr m 


( lr n i Q ~ ^ \ 


( lr m 

1 .EVlll O ( 


( lr m 




1987 Jun 09 


46955.7364 


-41.07 


1.75 


31.63 


4.57 


-4.71 


0.6551 


1988 May 29 


47310.7294 


-35.85 


0.81 


30.64 


2.13 


-6.27 


0.6305 


1988 Jul 03 


47345.7140 


-33.98 


1.24 


32.58 


3.23 


-5.53 


0.6976 


1989 Jan 19 


47546.0505 


-31.65 


1.75 


27.24 


4.57 


-3.74 


0.5348 


1989 Jan 23 


47550.0308 


-27.58 


1.94 


18.37 


5.08 


-4.11 


0.7700 


1989 Apr 12 


47628.9749 


-20.97 


1.87 


10.95 


4.89 


-11.84 


0.4345 


1989 Apr 16 


47632.9116 


-36.85 


1.40 


34.40 


3.66 


-6.95 


0.6671 


1989 Apr 19 


47635.9830 


-9.80 


2.47 








0.8486 


1989 May 12 


47658.9375 


19.60 


2.85 


-35.74 


7.47 


-5.70 


0.2049 


1989 May 15 


47661.8636 


-12.68 


1.81 


4.99 


4.72 


-0.68 


0.3778 


1989 May 18 


47664.8394 


-32.58 


1.46 








0.5536 


1989 May 22 


47668.8128 


-24.29 


1.49 


13.80 


3.90 


-1.60 


0.7884 


1989 May 23 


47669.8250 


-13.61 


1.17 








0.8482 


1989 May 24 


47670.7618 


2.43 


1.28 


-11.98 


3.34 


-4.65 


0.9036 


1989 May 25 


47671.7267 


12.96 


1.20 


-34.44 


3.14 


-7.22 


0.9606 


1989 May 26 


47672.7810 


25.38 


1.94 


-46.11 


5.08 


-5.68 


0.0229 


1989 Jun 25 


47702.7097 


-21.89 


1.81 


13.28 


4.72 


-3.30 


0.7913 


1989 Jun 26 


47703.6753 


-14.55 


2.64 








0.8483 


1989 Jul 12 


47719.6712 


-24.27 


1.87 


24.79 


4.89 


-4.34 


0.7934 


1990 Feb 06 


47929.0414 


24.93 


1.17 


-48.96 


3.05 


-7.29 


0.1644 


1990 Jun 03 


48045.7813 


28.00 


1.08 


-46.02 


2.82 


-1.28 


0.0621 


1990 Jun 11 


48053.7827 


-32.20 


1.20 


23.17 


3.14 


-2.70 


0.5349 


1990 Jun 12 


48054.7620 


-33.88 


1.22 


37.06 


3.19 


-9.89 


0.5928 


1991 Mar 31 


48346.9125 


-11.90 


1.35 


1.82 


3.52 


-11.89 


0.8550 


1991 Apr 01 


48348.0054 


4.70 


1.32 


-11.27 


3.46 


-1.14 


0.9195 


1991 Apr 02 


48348.9374 


18.10 


1.30 


-34.79 


3.40 


-7.49 


0.9746 


1991 Apr 25 


48371.9536 


-4.57 


2.85 


-4.59 


7.47 


-4.57 


0.3346 


1991 May 01 


48377.9066 


-36.38 


1.87 


41.07 


4.89 


-8.52 


0.6863 


1991 May 03 


48379.8833 


-20.94 


1.49 


14.65 


3.90 


2.12 


0.8031 


1991 May 21 


48397.8579 


-10.87 


2.85 


-1.15 


7.47 


-1.21 


0.8651 



Table 2 — Continued 



UT Date 


HJD 




0vl 


V2 




v tert 


SB 




(2,400,000+) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


Phase 


1991 May 23 


48399.7728 


15.13 


1.43 


-40.74 


3.74 


-2.29 


0.9783 


1991 Jun 03 


48410.8182 


-38.65 


1.24 


34.68 


3.23 


-5.73 


0.6309 


1992 May 12 


48754.8943 


12.72 


1.56 


-31.59 


4.09 


-6.36 


0.9612 


1992 Jun 09 


48782.7909 


-36.16 


1.24 


30.45 


3.23 


-6.97 


0.6095 


1992 Sep 06 


48871.6234 


-14.20 


1.65 


1.00 


4.31 


-5.29 


0.8583 


1993 Mar 31 


49077.9395 


29.73 


1.65 


-54.00 


4.31 


-8.79 


0.0488 


1993 Apr 03 


49080.9905 


14.98 


1.65 


-35.60 


4.31 


-5.52 


0.2290 



Table 3. Infrared Spectroscopic Observations and Radial Velocities 



UT Date 


HJD 






V2 


Cu2 


v tert 


SB 




(2,400,000+) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


(km s- 1 ) 


(km s _1 ) 


Phase 


2000 Jun ll a 


51706.8238 


-12.69 


2.54 


-4.26 


2.67 


-4.26 


0.3805 


2002 Jul 18 


52473.7874 


-38.14 


1.19 


25.86 


1.36 


-4.26 


0.6977 


2003 Feb 08 


52679.0828 


-20.60 


2.63 


4.90 


3.00 




0.8279 


2003 Apr 20 a 


52750.1016 


21.33 


1.79 


-51.03 


1.78 


-4.27 


0.0242 


2004 Jan 27 


53032.1879 


-38.17 


1.64 


25.83 


3.28 


-5.17 


0.6917 


2004 May 24 b 


53149.8809 


-38.32 


1.69 


29.78 


1.61 




0.6457 


2007 Apr 30 


54220.9531 


3.17 


2.10 


-28.82 


4.21 




0.9317 



a NIRSPEC+AO Observation; SB and tertiary spectra separated. 
b NIRSPEC+AO Observation; SB only. 



Table 4. SB Orbital Elements and Derived Properties 





Single-Lined 




Double-Lined 




Element /Property 


Visible Light a 


Visible Light 


Infrared 


Combined 


P (days) 


16.925 


16.9257 ±0.0010 


16.9231 ±0.0023 


16.9243 ± 0.0002 


7 (kms -1 ) 


-5.0 


-6.28 ± 0.25 


-9.99 ± 0.38 


-8.06 ± 0.29 


d7/dt (kms^ 1 yr _1 ) 


— 


-0.01 ± 0.18 


-0.04 ± 0.22 


-0.25 ± 0.04 


Ki (kms- 1 ) 


— 


33.06 ±0.38 


32.39 ±2.01 


33.12 ±0.33 


K 2 (kins" 1 ) 




42.70 ±0.89 


42.74 ±2.55 


42.43 ±0.67 


c 


0.10 


0.115 ±0.010 


0.120 ±0.018 


0.113 ±0.008 


uj (degrees) 




317.0 ±5.8 


322.1 ±33.6 


319.2 ±5.0 


T (MJD) 




52665.36 ±0.33 


52834.46 ± 1.60 


52698.92 ±0.23 


Mi sin 3 i 




0.421 ±0.022 


0.414 ±0.076 


0.417 ±0.017 


M 2 sin 3 i 




0.326 ±0.014 


0.314 ±0.059 


0.325 ±0.011 


/H b 


0.051 








q = M 2 /M 1 


0.71 ±0.05 c 


0.77 ±0.02 


0.76 ±0.07 


0.78 ±0.01 


a\ sin i (Gm) 


7.18 


7.644 ± 0.088 


7.482 ± 0.463 


7.659 ±0.077 


a 2 sin i (Gm) 




9.872 ±0.207 


9.875 ±0.589 


9.811 ±0.154 






0.969 


0.367 


0.891 


iVi ,N 2 


27 c , - 


37 , 33 


7,7 


44 , 40 


a vl , W v2 (kms -1 ) 




1.65 , 4.21 


1.94 , 2.56 


1.69 , 3.92 



Mathieul (119941 ) 



3 Mass function, f(m) = (m 2 sim) 3 / (m 1 +m 2 J 



iWalter et al.l (119941 ) identifies the three components and estimates the mass ratio from 27 
visible light observations. 



Table 5. H-R Diagram Derived Properties of NTTS 155808-2219 



Property Value 



SpT 

A,, (mag) 
log(Li/L ) 
log(L 2 /L ) 
log(L 3 /L ) 
M l (M ) 
M 2 (M ) 
M 3 (M ) 
M 2 /Mi 
Ms/Mi 

M t /Msb 
Ri (R ) 
Age (Myr) 



M2/M3/M2.5 
0.5 ± 0.2 
-0.56 ± 0.04 
-0.84 ± 0.04 
-0.66 ± 0.04 
0.50 ± 0.10 
0.38 ± 0.10 
0.44 ± 0.10 
0.76 ± 0.14 
0.88 ± 0.14 
0.50 ± 0.17 

1 qO+0.19 

o o +2.2 
00 -1.3 



Table 6. Tertiary Astrometric Measurements 



Date 


Position Angle 


Op A 


Separation 


@ sep 


(UT) 


(degrees) 


(degrees) 


(") 


(") 


1994.332 a 


313.7 


1.2 


0.193 


0.005 


1998.348 b 


317.5 


0.8 


0.21 


0.04 


2000.444 


315 


4 


0.22 


0.04 


2003.300 


321 


4 


0.23 


0.04 



Kohler et all fj2000h 



Woitas et all ( 120011 ) 
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Fig. I.— NIRSPEC H-band spectra of NTTS 155808-2219. The three epochs of AO 
observations are SB only, while non-AO spectra include the SB+Tertiary. All spectra have 
been corrected for heliocentric motion, flattened, normalized to a continuum level of unity, 
and separated in the y-axis by an additive constant. 



'. T eff =4209K Synthetic Template 




1.550 1.555 1.560 1.565 
Wavelength (/xm) 

Fig. 2. — Spectra of the NTTS 155808—2219 tertiary and cross-correlation templates. Top: 
Best matching synthetic spectrum calculated with the PHOENIX model atmosphere code 
and rotationally broadened to 25 kms^ 1 . Middle: Best matching observed template, GL 
752A, rotationally broadened to 18 kms^ 1 . Bottom: Averaged spectrum from the two 
NIRSPEC+AO epochs (UT 2000 June 11 and 2003 April 20) which isolated the tertiary. All 
spectra have been corrected for heliocentric motion, flattened, normalized to a continuum 
level of unity, and separated in the y-axis by an additive constant. Our analysis shows that 
the GL752A spectrum provides the highest correlation, while the synthetic template gives 
more consistent radial velocities and correlation values within 5% of the GL 752A results. 
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Fig. 3 . — Th e primary vs. secondary radial velocities for NTTS 155808—2219 following 



Wilson! ( 1194 ll ). The dotted line shows the weighted fit to the visible light data (bare error 



bars) and the dashed line shows the fit to the IR velocities (boxes). The mass ratio, q, 
is the negative of the slope of the fit and the center-of-mass velocity is determined by the 
equation 7 = (y-intercept)/(l + q). The vertical offset of the two fits gives a center-of-mass 
acceleration of d7/dt = -0.30 ± 0.07 (kms _1 yr _1 ). 



t 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




j i i I i i i I i i i I i i i I i L 



0.0 0.2 0.4 0.6 0.8 1.0 

Phase 

Fig. 4. — Radial velocity minus the center-of-mass velocity as a function of phase for 
NTTS 155808—2219. The diamonds represent the primary star data, and the triangles the 
secondary star data. Filled symbols are IR velocities and open symbols are visible light ve- 
locities. The best fit orbital solution to the combined visible light and IR data is represented 
by a solid line for the primary star and a dashed line for the secondary star. 
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Fig. 5. — Change in the center-of mass velocity of the binary as a function of time. Each point 
represents the weighted average of the primary and secondary residuals after subtracting the 
orbital motion. The dashed line corresponds to the fitted 7 velocity from our orbital fit as a 
function of time. The center-of-mass acceleration from the orbital solution is d7/dt = -0.25 
± 0.04 (kms _1 yr _1 ), which is consistent with the method in Figure [3] 




Fig. 6. — T e jj and luminosity for the NTTS 15580 8 —221 9 components plotted on the pre- 
main sequence evolutionary tracks of iBaraffe et al.l (119981 ) for solar metallicity and mixing 
lengths a = 1.0 for M < O.5M and a = 1.9 otherwise. Luminosities for the three compo- 
nents were derived from the total 2MASS magnitudes for the s ystem and the s pectr oscopic 
component flux ratios. Effective temperatures were chosen from lLuhman et al.l (120031 ) for an 
M2 ± 1 primary, M3 ± 1 secondary and M2.5 ± 1 tertiary. The mass tracks and isochrones 
are labeled. For objects lying between labelled isochrones and mass tracks, the ages and 
masses were estimated by interpolation. The three NTTS 155808—2219 components appear 
to be coeval with an age of 3.3 tU Myr. 
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Fig. 7. — Masses of the NTTS 155808—2219 SB components as a function of orbital incli- 
nation. The values of Misin 3 i and M2sin 3 i derived from the orbital solution are shown with 
la uncertainties (c urved dashed and das h-dot lines). Primary and secondary mass estimates 
from the tracks of IBaraffe et al.l (119981 ) are over-plotted (horizontal dashed and dash-dot 
lines), also with la uncertainties. The inclination must be > 61° for the model masses to be 
consistent with the minimum masses within the la uncertainties. 
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Fig. 8. — Radius of the NTTS 155808—2219 primary as a function of stellar inclination. 
The primary star's radius, derived from the luminosity and temperature used in Figure [6j 
is over-plotted with la unce rtainties (cros s -hatch ed region). We determine Rsim from the 
rotational period reported by lAdams et al.l ( 119981 ). 4.3 ± 0.3 days, and our measured vsini's. 
Our IR spectra show the highest correlation with templates rotated to vsim w 20 kms" 1 ; our 
visible light data yield vsim < 10 kms^ 1 . Our radius derived for the primary is consistent 
with values of Rsim for a range of vsim of 10—20 km/s if the stellar inclination is >33 
degrees. 



